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Abstract. This paper is devoted to the problem of sample path large deviations for multidi- 
mensional queueing models with feedback. We derive a new version of the contraction principle 
where the continuous map is not well-defined on the whole space: we give conditions under 
which it allows to identify the rate function. We illustrate our technique by deriving a large 
deviation principle for a class of networks that contains the classical Jackson networks. 
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Introduction 



This paper is concerned with the theory of large deviations of stochastic processes related 
to discrete event systems. As opposed to classical stochastic dynamical systems, for which 
the evolution is continuous and described by a stochastic differential equation, discrete event 
systems are characterized by synchronization mechanisms that prevent most of the classical tools 
to apply. We present here a new approach for the analysis of the sample path large deviations 
of such processes. Unlike standard methods that require establishing upper and lower bounds, 
our method relies on the analogy between the theory of weak convergence and the theory of 
large deviations. This analogy is well-known and has been studied by many authors, we refer 
to the recent book of Feng and Kurtz jllj that surveys this field. We should in particular 
quote the work of Puhalskii 20 , [221 quite similar to our approach. However, we will not use 
the framework of idempotent measures developed by Puhalskii. We discuss in more details our 
general methodology and its relation with the existing literature after the description of the 
queueing networks we consider. 

To apply our method we choose a class of queueing networks with Bernoulli routing, where 
feedback is allowed. The discontinuous dynamic of queueing networks makes it hard to study 
and large deviations results in the literature are treated on a case by case basis as in the work of 
Ganesh and Anantharam |T^, Bertsimas, Paschalidis and Tsitsiklis 131 or Ramanan and Dupuis 
j24j . As we will see, adding the possibility of feedback makes the problem much harder. For 
queueing networks with feedback, existing large deviations results are restricted to networks 
described by finite-dimensional Markov processes, see the works of Dupuis, Ellis and Weiss jlOj . 
Dupuis and Ellis 0| and Ignatiouk- Robert jl5| . In this paper, we consider networks where 
the output process of a queue is modeled by a reflection mapping. This class contains the 
classical Jackson networks and our large deviations results extend existing results for this class 
obtained by Atar and Dupuis |2j and Ignatiouk- Robert I14J. Our technique allows to obtain 
large deviations results under non-exponential assumptions. This case corresponds to networks 
with autonomous service and gives an approximation for queueing networks where each station 
acts as a standard single server queue. While preparing this paper, the author became aware 
of the work of Puhalskii |^ who considers generalized Jackson networks. The form of the rate 
function for the queue length process obtained in coincides with our result, which confirms 
the intuition that in the large deviation regime, networks with autonomous service approximate 
well generalized Jackson networks. We will discuss more carefully this result in Section [2.31 

In the next section, we give an overview of the general methodology and then introduce the 
general notation. Section ^ gives an extension of the contraction principle that will allow us to 
identify the rate function. Our result is stated without any reference to any specific discrete 
event system and could be applied to other systems. In Sections |31 and |1J we apply our method 
to the case of a queueing network. 

General methodology. For simplicity, we adopt here the notation corresponding to our ex- 
ample of queueing network. As in jl6j or ^7j, we define the arrival and departure processes A 
and D of each station of the network as the solution of the fixed point equation 



Here Net is a process that describes all the primitives of the network as service times at the 
different stations, routing decisions, arrival times in the network. The maps F and $ describe 
the dynamic of the network. 



(0.1) 
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We consider a sequence of queueing networks {Net„}„, where the primitives are counting 
processes (i.e. Net„ belongs to a space denoted by f ) and satisfy a large deviation principle 
(LDP). We denote by l'^''* the corresponding rate function. It is known that the map ^' is well 
defined if the primitives of the network are counting processes, see j7j or ^Hj and we denote 
(A„,D„) = ^'(Net„). It is natural to ask whether ^ is well defined for processes in D, the 
space of cadlag non-decresing functions or at least for absolutely continuous functions. If this 
was true, and if ^ was shown to be continuous, then we would get thanks to the contraction 
principle that the sequence of processes {(A„,D„)}„ satisfies a LDP with good rate function 

(0.2) "l^'°(A,D) = inf {l^'=*(Net), ^(Net) = (A,D)}." 

However, the map ^ turns out not to be well defined for all possible limits of a sequence of 
networks {Net„}„ G £^ as defined previously. In particular, the fixed point equation (|U.1() can 
very well be stated for processes in D but then may have several different solutions as noted by 
Majewski We give in the appendix a simple example. 

To circumvent this difficulty, we adopt the following strategy. We find a domain PNet C £ 
satisfying the following constraints: 

• the map ^ is well defined on I^Net; 

• any solution (A, D) of the fixed point equation (jU.lj) associated with a "continuous" 
Jackson network Net can be approximated by a sequence {Net„} G P^g^ such that 

(0.3) Net„ Net, 

(0.4) ^'(Net„) ^ (A,D), 

(0.5) l^'=*(Net„) ^ |N''*(Net). 

Hence in order to remove the quote from (|U.2() . we follow a quite standard method of proofs 
for large deviations of stochastic processes analogue with the theory of weak convergence 
it consists of first verifying a compactness condition and then showing that there is only one 
possible limit. In our context, we proceed as follows: 

(1) we show that our sequence of processes is exponentially tight; 

(2) we use PNet to determine the rate function. 

In Section Q we give the theoretical framework that shows how any domain verifying assump- 
tions (jnSl), and (1031) det ermines the rate function. This result is stated in great generality 
(without any reference to our specific problem) and could be of independent interest since this 
method of proof could be applied to other dynamical systems (with discontinuous statistics) . 

Notation. For {E,d,<) a complete, separable metric space with partial order <, we denote 
by D(i?) the space of cadlag non-decreasing i?-valued functions defined on M_|_ with Skorohod 
(Ji) topology and by C{E) the space of continuous non-decreasing E'-valued functions defined 
on M+. Restricted to C{E) the Skorohod topology is just the compact uniform topology. 

For x,y (z M'^^, we write x < y \i x^*-* < y*-*-* for all i. We denote by A the minimum and 
by V the maximum in M.^ . For X, Y G B(]R^), we write X < Y if X(t) < Y(t) for ah t > 
and for maps F, G G D(Mf)2, we denote F < G if F(X) < G(X) for all X G B(]R^). For 
X G M^, we denote = vf^^x^ and for X G D(Mf ), we denote ||X|| = sup^ ||X(f)||. We 
denote Bo(F) = {/ G B(F), /(O) = 0} and Co(F) = {/ G C(F), /(O) = 0}. 

A piecewise linear function is a continuous function such that there exists a partition r = 
(to = < ti < . . . ) with tk ^ oo and such that the function is linear on each interval (t^, tfc+i). 
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For any function / G D(M^'), we define the polygonal approximation of / with step 1/n as the 
(piecewise linear) function 



n I \ \ n I \ n 



is the set of substochastic matrices of size K x K. For M G M-'^, we denote by p{M) 
its spectral radius, by M* its transpose and M^*) denotes the line M^*) = (M(*'^\ . . . M^*'^)). 
In particular, we will identify a function P € B(M-^) with its K components P^*^ S D(M:^), 
where pW(t) = (p(^'i)(t), . . . p(^'^)(t)) with Y.j^^''^\i) < 1 for ah t > and all i. Note that 
for M, iV G M^, we have M < if M^*'^) < iV^*'-') for all i and j. 

We will use the Kullback-Leibler information divergence, which is a nonsymmetric measure 
of distance between distributions in the sense that for any two distributions P and R on 
where is a finite set, 

D{P\\R)= P(x)log^^^^' 



R{x] 



is nonnegative and equals if and only if P = i?. We use the standard notational conventions 
logO = — oo, log ^ = oo and OlogO = Ologg = 0. For any fixed R, the divergence D(P||i?) is a 
continuous function of P restricted to {P, S[P) C S{Q)} where S{P) denotes the support of P 
(see m). 

For P G M^, we denote by P the K x [K + 1) stochastic matrix obtained as follows: for all 
i,j < K, P(^'J) = and = 1 - ^f^^ For P,Re M^, we wiU denote 

D{P\\R) := D(P||i?) 

= y log ( ^] +y(i-y p(^A log ( ^-^^^] 

K 



^D(P«||ii«). 
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1. An extension of the contraction principle 

Let E,T be complete separable metric spaces. Let G:fxJ^— s-Mbea continuous function. 
We assume that there exists 2? C such that for all x G there exists an unique y ^ T such 
that y) = 0. We denote it by, y = H{x) where H -.V J^, 

Vx G V, G(x, y) = ^ y = H{x). 

Proposition 1.1. Let {Xn}n be a sequence of £- valued random variables and {Yn}n be a se- 
quence of !F -valued random variables. We assume that each sequence is exponentially tight. As- 
sume that the sequence satisfies a LDP with good rate function and that G{Xn, Yn) = 
a.s. for all n. 

We assume that for all {x,y) such that G{x,y) = and \-^{x) < oo, there exists a sequence 
Xn X, such that x„ G P for all n, H{xn) — > y and \^ {xn) \'^{x). We denote by S{x,y) = 
{xn}n this sequence. IfG{x,y) ^ or \^ (x) = oo, we take S{x,y) = and we denote S{y) = 
Ux{S{x,y)} (which might be empty). 
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Then the sequence {X„, Y"„}„ satisfies a LDP with good rate function: 

(1.1) |^.^(x,y):=| G[x,y) = 0, 

[ cxD, otherwise. 

In particular, if € V for all n and if the sequence {H{Xn)}n is exponentially tight, then 
it satisfies a LDP in T with good rate function: 

(1.2) l^^'')(y) := mf{ lim l^(x„), E S{y)}. 

n— >oo 

Remark 1.1. • There are alternative ways of expressing the rate function, 

l™(y)=inf{l^(x), ye//-}, 

where := {y G J^, 3x„ — > x, H{xn) y}. l^^-^^ is the lower semicontinuous 
regularization of the following function defined for y G H{D) C T, 

FW(y) :=inf{l^(x), y = /7(x)}. 

The main interest of the definition (|1.2j) is that the rate function is computed only thanks 
to the sequences S{x,y) G V^. 
• Note that if LL{D) is closed (in particular if P = £) then this proposition follows from 
the contraction principle (for an extensive discussion of this principle, see the work of 
Garcia ^Hj)- Roughly speaking. Proposition 11.11 tells us that if D is dense in a certain 
sense in £, then the contraction principle still holds for the map H. 

Proof. Thanks to Lemma 3.6 of 11 , the sequence {X„,y„}„ is exponentially tight. Then by 
Theorem 3.7 of ^T], there exists a subsequence {n^} along which the sequence {X„^,y„^}„^ 
satisfies a LDP with a good rate function. If we can prove that there is a unique possible rate 
function (that does not depend on the subsequence {n^}) then the proposition will follow. 

Hence, for simplicity of notations, we still denote the extracted subsequence by {Xn,Yn}n 
and we assume that {Xn,Yn}n satisfies a LDP with good rate function 1^'^. We will show that 

Consider the continuous mappings Hi and H2 from £'x.7^to£'x.Fx]R, 

Hi{x,y) := {x,y,G{x,y)), H2{x,y) := (x,y,0). 

We have clearly Hi{Xn,Yn) = H2{Xn,Yn) a.s. Moreover thanks to the contraction principle, 
{Hi{Xn,Yn)}n and {H2{Xn,Yn)}n Satisfy LDPs with the good rate functions 

1^^ (x, y, z) = inf{F'^(x, y),z = G{x, y)} F^x, y, z) = inf{F'^(x, y),z = 0}, 

where inf = 00. Since Hi{Xn, Yn) = H2{Xn, Yn), we have Fi = F2. Now we have, 

F^(x, y) = inf{Fi (x, y, z)} = inf{F^(x, y), G(x, y) = 0}, 

hence F'^(x,y) = 00 as soon as G{x,y) ^ 0. It remains to show that G{x,y) = implies 
F'^(x,y) = F(x). We have clearly F(x) < F'^(x,y) for all (x,y) since satisfies a LDP 

with good rate function 

F(x) = inf{F^(x, y), y G ^, G(x, y) = 0}. 
In particular, the definition of D implies F(x) = F'^(x, H{x)) for x G P. 
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Take {x,y) such that G{x,y) = and < oo. There exists x* ^ x with x* G V, 

i?(x*) — > y and ^ l^(x). Thanks to the lower semicontinuity property of 1"'*-'^, we can 

find for any (5 > 0, an e > such that 

^A(F'^(x,y)-5) < inf F'^(x,z), 

where B{y,e) is the closed ball in of center y and radius e. 

Thanks to the lower semicontinuity of the function x i^izeB{y,e) \^'^ {x, z), we have 

inf ~\^'^{x,z) < liminf inf F''^(x„,z) 

z£B{y,€) Xn^x zeB(y,e) 

< liminf inf V'^{x*„,z) 

n^oo zeB{y,e) 

< hm F(x:) = F(x), 

because H{x'^) S B{y,e) for sufficiently large n. Hence we proved that for any 5 > 0, | A 
(]'^'^{x,y) — 5^ < F(x) for {x,y) such that G{x,y) = and F(x) < oo, this concludes the 
proof of (HHJ). 

The various expressions of are now quite easy to obtain from 

(1.3) FW(y)=inf{F(x), G{x,y)=0}. 

For H1.2|) . note that since the set {x, G{x,y) = 0} is closed the minimum in (if it is finite) 
is attained for a certain x* with G{x*,y) = and F(x*) < oo. 

We prove now that 

inf{F(x), yeH'^} = inf{F(x), G{x,y) = 0}. 

If y G H^, then there exists x„ ^ x such that H{xn) — ^ y- Hence by continuity of G, we have 
G{x,y) = 0. Now if G{x,y) = and F(x) < oo, it follows from the assumptions that y G H^. 

To see that the last expression in Remark II. II is true, we show that for any open set O C 
we have, 

(1.4) inf F(^)(y) = inf {F(x), y = H{x)}. 

y&O yeO 

For y £ O and any x such that G{x, y) = 0, there exists Xn x, such that H{xn) y and 
F(x„) —>■ F(x). Hence for n sufficiently large, we have H{xn) G O and then 

inf {F(x), y = H{x)] < inf F(x„) < F(x). 

3/GO n 

Taking the minimum over all x such that G{x, y) = gives the > inequality in p.4() . the converse 
inequality is obvious. □ 



2. QUEUEING NETWORKS WITH BERNOULLI ROUTING: DESCRIPTION AND LARGE DEVIATIONS 

RESULTS 

2.1. General setting and notation. We start with the basic model for an isolated queue and 
refer to ^ for more details on the relationship with other models of the literature. 
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The model for an isolated queue is in term of two primitive quantities belonging to ©(M^): 
the arrival process A and the service process S. The departure process D is a derived quantity 
that is obtained as a functional of the arrival and service processes as follows: 



(2.1) D(t) := inf {S(t) - S{s) + A{s)} A S(t). 

0<s<t 



From a mathematical point of view, if 7^ : D ^ D (where B is the space of cadlag M- 
valued functions defined on M-|_) is the one-dimensional Skorohod's reflection map defined by 
7^(X)(^) := supo<3<i {X(t) - X(s)} V X{t). We have D = A - 7^(A - S). It is easy to see that 
D € B(M+) and D < A. 

The queue length process is defined as the difference of the arrival process and the departure 
process, 



q{t) := A{t) - D(t) = sup {A(t) - A(s) - (S(t) - S(s))} V (A(t) - S(t)) . 

0<s<t 



If the arrival process A and the service process S are counting processes, this model is called 
a single queue with autonomous service: the queue length is increased by one whenever there is 
an arrival from the arrival process and the queue length is decreased by one whenever there is 
an arrival from the service process and the queue is not empty (see ) . Note in particular that 
in the case where the process S is a Poisson point process, then this model is a standard -/M/l 
queue. 

We now consider networks obtained by interconnecting queues modeled by ()2.1() when the 
departure process of one queue is randomly routed to the other queues as for Jackson networks. 
The networks we consider are characterized by the fact that service times and routing decisions 
are associated with stations and not with customers. This means that we associate to each of the 
K stations three predefined counting processes: an arrival process, a service process and a routing 
process. The arrival process and the service process of station k are described by the sequences 
of exogenous arrival times {7j'^^}j>i and service times {cj'^^}j>i- If there is no exogenous arrival 

(k) 

at station k, we use the convention Tj = oo for all j. When the j-th customer has completed 

(k) / (k) 

his service at station k, he is sent to station Vj (or leaves the network i/^ = K + 1) and 
is put at the end of the queue on this station, where {z^j'^''}j>i is also a predefined sequence, 

called the routing sequence. The sequences {Tj''^}j>i, Wj''^}j>i and {i^j''^}j>i, where k ranges 
over the set of stations, are called the driving sequences of the network. A network will be 
defined by {{trf ^}j>i, {z^f ^}i>i, {Tj''^}j>i, n^^\ 1 < A; < if}, where {n^^\ . . . describes 
the initial condition. The interpretation is as follows: at time t = 0, in node /c, there are 
n'-'^) customers with service times , • • • , o"^^!:) (if appropriate, a^"^ may be interpreted as a 
residual service time). In particular at time 0, the total number of customers in the network is 
^(1-^) =n« + ...nW. 

In what follows, we will describe the driving sequences thanks to their associated counting 
functions. We will use the following notation: (7*^'^)(l,n) = X]j=i ^ f°^' Q < k < K . 
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We define the sequence of networks Net„ = {S„(t), Pn(t), N„(t)} with 



N«(i) = 


\ k 


s®(0 = 


1 ^ 

-Z^l{<TW(l,fc)<nO 




k<nt 



Note that we allow the initial queue length to depend on n, N,*'' (0) = n^^ but the other driving 
sequences describing the arrival times, the service times and the routing decisions do not depend 
on n. Note also that if there is no exogenous arrival at station i, we have N^n\t) = Ni'^(O) for 
ah t > 0. 

For the network Net„, we denote the corresponding input and output processes of each 

(k) (k) 

queue k of the network by A„ and respectively. We will use the following notation 

An = {An^ , . . . , A^^^ ) and D„ = (D^^'' , . . . , D^^-* ) . We now describe how the processes A„ and 
D„ are obtained form Net„. 

We define the map F : Do(K+) x Bo(M'^) x D(Mf ) D(Mf ) as follows: 

K 

r(X, P, N) W (t) := N« {t) + P(^'^) (X(^') {t)). 

The following lemma is straightforward. 

Lemma 2.1. The map F is continuous for the compact uniform topology and non- decreasing in 
its first argument. 

We define the map $ : D(Rf ) x Do(Mf ) ^ Bo(Mf ) as follows: 

$(X,Y)«(i):= inf |Y«(t)-Y«(s)+X«(s)) AY«(t). 
o<s<t L J 

Lemma 2.2. The map $ is continuous for the compact uniform topology and non- decreasing in 
its first argument. 

Proof. We can clearly consider the map $ with K = 1 only. Let TZ be the one-dimensional 
reflection map, we have $(X, Y) = X — TZ{X. — Y). It is easy to see that for any T > 0, 

sup |7^(x)(^) - 7^(x')(^)| < 2 sup \x{t) - x'{t)\, 

0<t<T 0<t<T 

from which the continuity of $ follows. Its monotonicity is obvious. □ 

Remark 2.1. Consider the mapping $ with K = 1 and Y(i) = ij,t, with > 0. If // = 0, 
since $(X,Y) < Y, we have $(X,Y)(t) = for ah t. If /x 7^ 0, we have $(X,Y)(t) = 
info<s<t{X(s) + n{t — s)}. Moreover if X is a concave function, then this equation reduces to 
$(X, Y)(i) = X(t) A fit. Hence we can write 

Y{t) = lit, with II >0 => $(X, Y)(t) = litA inf {X(s) + n{t - s)}, 

0<s<t 

if moreover X is a concave function $(X, Y)(t) = iit A X(t). 
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(2.2) 



It is easy to adapt the proof of Theorem 2.1 of or Proposition 2.1 of to show that the 
following fixed-point equation: 

A„ = r(D„,P„,N„) =r(D„,Net„), 
D„ = $(A„,S„) = $(A„,Net„), 

has an unique solution when each component of nS„, nP„ and nN„ is a counting function (i.e. 
non-decreasing function of D(M^) or ©(M-'^) that is piece-wise constant with jumps of size one). 
In this case the corresponding functions nA„ and nD„ are also counting functions and we denote 
the solution of by ^'(S„,P„,N„) = ^'(Net„). 

Remark 2.2. Note that the only difference between our model and generalized Jackson net- 
works as described in JB] resides in the queueing mechanism (|2.1|) which is sometimes called 
autonomous. Consider a network Net = {S,P,N} where the processes are counting processes. 
Then due to some monotonicity arguments, it is possible to relate (see jH]): 

- the process (A,D) associated to Net with the dynamic described in |16j : 

- the processes ^'(Net) = (A,D) solution of the fixed point equation. 

Note that in the case where the process S is a Poisson point process, our model is exactly a 
Jackson network (see j2]). 

2.2. Stochastic assumptions. In what follows, it will be important to distinguish the nodes 
of the network that do not receive any exogenous customer, i.e. the nodes i €z S'^ with S = 
{i, T^'^ < oo}. A network Net = {S,P,N} is an object in £" C ©0(1^^) x Bo(M^) x inim^^ 
with the additional constraints: 

(1) NW(t) = N«(0) for all t, for i ^ S; 

(2) for ah < u < n, we have Ylf=i P'^''^Hu) - P^''^^v) <{u-v). 
Note that £ is closed in Bo{Rf) 

We define for (s^^), . . . , s^^)) G and (n^^), . . . , n(^)) G , the functions 



|S(s«,...,sW) = X]l^*''(5«), 

1=1 

l^(n«,...,n(^)) = El''"(-^^^) + -l{n«>o.^5}, 

ieS 

where each 1^^'^ (resp. l'^''' for i G 5) is a [0, cx3]-valued convex good rate function, attaining zero 
on M+ admitting a unique minimum at the point (resp. A*^*^ for i G 5) and with a domain 
open on the right. 

We assume that the sequence Net„ = {S„(t), P„(t), N„(t)} satisfies a LDP in the space £ 
with a good rate function l'^'^* given by 

/•OD 

(2.3) |N^*(S,P,N) :=|0(N(0))+ / l^(S(t)) + D(P(t)p) + |N(N(t))dt, 

Jo 

if the argument functions are absolutely continuous and equal to infinity otherwise. 
We make the following assumptions on the matrix R: 

(1) We assume that p{R) < 1. 
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(2) We assume that for all 1 < i < -ftT, we have 

(2.4) (A/' + A/'i2+ ••• +AAii^)» > 0, 
where M is the line vector of defined by AA^*) = Ijjg^}. 

We show now that our stochastic assumptions cover the case where S^'^) and N(*) (with i G S) 
are independent and correspond to renewal processes and where the routing is a Bernouilli 
routing associated with the matrix R that satisfies previous assumption. 

We recall here some results of Puhalskii [25 concerning large deviations of renewal processes 
and show that our assumptions on the rate function (|2.l'{|) are satisfied in the i.i.d case. Denote 
by {Ci) * ^ 1} 3- sequence of non-negative i.i.d. random variables with positive mean. Let 

a{e) = logEfe^^i 

e* = sup{0 > 0, a{e) < oo}, 

(2.5) a*{x) = sup{^x-a(e)} = sup{0x-a(0)}, 

e e<e* 

g{x) = xa*(l/x) = sup {0 — xa(^)}. 

6»<6»* 

Note that the function a is a convex function and differentiable on (— oo,0*) with a'(0) = 
E[(^i] > 0. In particular, we have Ivchq^q* a{9) = oo, from which we get the equality in (|2.5j) . The 
functions a* and g are convex rate functions. Introduce the sequence of processes {Cn}n- 

i 

Then Theorem 3.1 of 21 gives: If P(Ci > 0) = 1, then the sequence {C„}„ satisfies a LDP in 
B(M_i_) with the good rate function 

|C/ X _ f g{^{t))dt, if x G C(IR+) is absolutely continuous, 
^ ' 1 oo, otherwise. 

It then follows that g is a good rate function. Moreover, we have essinf = if and only if 
g(x) is finite for all x > E[Ci]~^ (note in particular, that in this case, the domain of g is open on 
the right). The proof of this fact can be found in |19j and follows the argument: from a{9) > 
E[Ci]e, we have g (E[Ci]"^) = and for all x > E[Ci]-\ we have g{x) = supg<Q {9 - xa{e)}. If 
essinf = 0, we have for arbitrary e > and for 9 <0, 

a{9) = logE le^'^A >ee + logP(Ci < e), 

hence for x > e~^, we have g{x) < — logP(Ci < e)- It is clear that if essinf d > 0, then for any 
X > essinf C-f"^, we have g{x) = oo. 

Concerning the large deviations of the routing processes given in term of the Kullback-Leibler 
information divergence, it follows directly from Corollary 6.1 of 1231 in the case of Bernouilli 
routing, i.e. when the sequences {i^j''^}j>i are sequences of i.i.d. random variables in [1,K] and 
independent in k such that 

P(z.f) = i) = i?(M. 
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2.3. Sample path large deviations for the queue length process. We now return to the 
sequence of queueing networks defined in Section |2j Recall that (A„,D„) correspond to the 
arrival and departure processes from each station. We now give our theorem for the queue 
length process defined as Qn = A„ — D„. 

Theorem 2.1. The sequence of processes {Qn}n satisfies a LDP in D(M:^) with good rate 
function that is finite for Q absolutely continuous given by: 

"°(Q(o)) + i§(o)(Q)^ 

where for q > 0, is a good rate function that is finite for absolutely continuous Q such that 
Q(0) = q and given by: 

POO 

|Q(Q):= / HQ(Q(s),Q(s))<is, 
Jo 

where is given by, 

HQ(Q,Q):=inf| J] 1^'' P«)l|z.(o>.(0} + E 1^" P^^^) + E + |N(iV) 

where E{Q) = {i, Q^*-* = 0} and the infimum is taken over the set of [D, P, iV) G Mf x x Mf 
such that 

Q = N + {P^ - Id)D. 



In jn], Puhalskii obtains a LDP for the queue length process of a generalized Jackson net- 
work with a rate function that coincides with Theorem 12.11 Note that our model is slightly 
different here since we model the dynamic of a queue by a reflection mapping. Still in the case 
of Poisson processes for the inputs, both models correspond to the (exponential distribution) 
Jackson network. Recall that the rate function for a Poisson process of rate A is given by (we 
keep the same notation as in l2.2() . 

(2.6) l'^(x)= r X±{t)log^ -±{t) + Xdt, 

Jo X 

for absolutely continuous functions x € C(M-|_). Hence if we replace (|2.6|) in the expression of I^, 
we obtain the rate function for the large deviations of a Jackson network. In this specific case, 
the rate function has been obtained in different forms by Atar and Dupuis 3 and Igniatiouk- 
Robert Jl] and some bounds have been computed by Majewski Compare to these results, 
our representation has the advantage of being quite intuitive, in the sense that each term is easy 
to interpret. If we interpret D, P, N as instantaneous departure, routing and exogenous arrival 
rates, then + (P* — Id)D is just the vector of rates at which the queue lengths vary. Hence 
given a rate of change of Q, the system behaves in such a way to minimize the instantaneous 
"costs" of departure, routing and exogenous arrival rates over all the rates that yield the desired 

Q 

From a methodological point of view, the argument of 19 is quite different from ours since 
the density condition (that we could compare to our Proposition II. 1|) is verified on the rate 
function 1^ (see condition (D) in jl^) whereas we are checking the density argument on the rate 
function of the inputs. 
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3. Extension of ^ to piece-wise linear networks 

In this section we consider processes that are continuous, i.e. in C{E), hence topological 
concepts refer to the compact uniform topology. 



We first recall Proposition 3.2 of [T 

Proposition 3.1. Given a K x K substochastic matrix P with p{P) < 1 and vectors {a,y) € 
M.^ , the fixed point equation 

has a unique solution x(y,P,a). Moreover, {y,a) x{y,P,a) is a continuous non- decreasing 
function. 

We first consider a linear network Net and show that the mapping ^ (defined as the solution 
of the fixed-point Equation (|2.2|l ) is well defined for such a network. By linear, we mean the 
following N(*)(t) = A^(^) + X^'h, with A^*) > and N^'^ G M+, S^'\t) = i^^'h, with ij.^'^ > 0, and 
p(»j)(t) = p(i'^)t. We assume that p{P) < 1. 

Lemma 3.1. Under previous assumptions, the fixed point equation has an unique solution 

Xf[p,P,N,X]{t) =x{nt,P,N + Xt), where fi = iV = (AtW). andA = (A«)i. 

Proof. Since /i, P, N, X are fixed here, we omit to explicitly write the dependence in these vari- 
ables. In this case, the fixed point equation (|U.1|) reduces to (see Remark l2.ip 



(3.1) 




A info<,<t{A«(s) + /i«(t - s)}. 



(3.2) 



Thanks to Proposition 13. H Xj(t) = x{pt, P,N + Xt) is the unique solution of the fixed point 
equation 

f A«(t) = + A«t + J2f=i P^^''W\t), 

I D«(t) = A«(t)A/i«t. 

We prove now that is the unique solution of the fixed point equation 1)3. 1|) . 

For simplicity, we denote the fixed point equation (|3.1|) . resp. 1)3. 2() . by A = F{A), resp. by 
A = F{A). Note that these functions are non-decreasing, continuous and such that F < F. 

From < X/, we get < F(0) < F(0) < F(X/). Hence F"(0) / L < X/ and F(L) = L. 
Moreover for any solution Y of the fixed point equation (|3.1() . we have L < Y < Xj because 
Y = F{Y) < F{Y) and F"(Y) / X/. 

Since is a concave function, we have F(0) = F(0) and hence it is still a concave function. 
Hence we have F^{0) = F"(0) since the image by F of a concave function is a concave function 
and F = F on the subspace of concave functions. Hence we have L = Xy which concludes the 
proof. □ 

In order, to extend ^' to piece-wise linear networks, we proceed step by step on each in- 
terval where the driving functions S,P,N are linear. The following lemma allows to glue the 
constructed solution on each adjacent interval. 
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Lemma 3.2. Let A, S G D(M+) x Bo{R+) and D = ^>(A, S). Define A, S G D(M+) x Bo(R+) 
as follows 

A{t) := A(t + 'u) -D(u), 
S{t) := S{t + u)-S{u). 

Let D = ^{A, S), then we have 

D(t) =D(t + ti) -D(n). 

Proof. We show that for D = $(A, S), we have 

D(t + n)-D(u)= inf {S(t + u) - S(s) + A(s) - D(u)} A {S(t + n) - S(u)} , 

u<s<t+u 

from which the lemma follows. 
We write 

D(t + n)-D(n) = inf {S(t + u) - S(s) + A(s) - D(n)} 

0<s<u 

A inf {S(t + u)-S(s) + A(s)-D(n)}A{S(t + u)-D(u)}, 

n<s<t+n 

Since D(ti) < S(n), we have to prove that 

S{t + u)-S{u)> inf {S(t + u)-S(s) + A(s)-D(n)} A{S(t + n)-D(u)}. 

0<s<n 

This will follow from, 

inf {S(t + n)-S(s) + A(s)-D('u)} = S{t + u) - S{u) + inf {S(u) - S(s) + A(s)} - D('u) 

0<s<u 0<s<u 

< S{t + u)-S{u). 

□ 

We consider now piece- wise linear networks: the functions u N(*)(ii), and 
u 1-^ p(*'-')(ti) are continuous piece-wise linear functions such that N'^*)(0) G M+ and S(*'*(0) = 
p(»j)(0) = and p{P{t)) < 1 for all t > 0. 

Proposition 3.1. For a piece-wise linear network, there exists an unique solution of the fixed 
point equation iP.l]) . We still denote by ^ the mapping that to any piece-wise linear network 
Net associates the corresponding couple (A, D). 

Proof. The existence is a direct consequence of monotonicity properties and continuity of the 
maps r and We define the sequence of processes {A [A;], D[A;]}fc>o with the recurrence equation: 

f A[A; + 1] = r(D[A:],Net), 
\ D[A: + 1] = $(A[A; + l],Net), 

and with initial condition D[0] =0. By the monotonicity properties of $ and F, we have 

< A[l] ^ ^>(0, Net) = = D[0] < $(A[1], Net) = D[l] 

r(D[0],Net) = A[l] < r(D[l],Net) = A[2], 

and the sequence {A[A;], D[A;]}fc>o is increasing. Note that D[A;] < S and hence the following 
limits are well defined 

lim A[k] = A and, lim D[A;] = D. 
Since T and $ are continuous, (A, D) is a solution of the fixed point equation (|0.1|l . 
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We now prove uniqueness. First recall that we call a, a partition of M-|-, any increasing 
sequence of points a = {a„}„ with = and o„ — > oo. For two partitions a = {an}n and 
/3 = we say that 7 = {gn\n is the union of a and /? if 7 is a partition such that for all n 

there exists m such that either = am or gn = bm- 

Let r = {tn}n be the union of the partitions associated with each function S, P, N. We define 
for X € M+, d{x, r) = min„{t„ — x, tn > x} > 0. 

Assume that we are given two solutions of the fixed point equation (|fl.l|) : (Ai,Di) and 
(A2,D2). First note that thanks to Lemmas 15.11 and I5.2[ any solution of (|n.lj) is absolutely 
continuous. Let z = inf{t, Ai(t) 7^ A2(t)}, in particular, we have Ai(t) = A2{t) and Di(t) = 
D2(t) for all t < z. 

Define u = mini d(D^\z) , t) A d(z,T) > 0, where the notation , can be replaced either by 1 
or by 2- We have that for t £ [0, u], 

S»(t) := S(*)(z + t)-S«(z) = t/i», 
p{^.i)(t) := p('j\-D^:\z) +t)- p(*'^')(Di*)(z)) = tP'-''^\ 

N» := n(*)(z + t)- N»(z) + a'^\z) - B^J\z) = tA^ + A'i\z) - B^J\z), 

Let A{t) = Xj[/^, P, A,(z) — 0,(2;), X]{t) be the unique solution associated to the infinite horizon 
linear network defined above. The associated departure process is D(t) = A(t) A fit. Let 
V = inf{t, infj t)^^\t) = u}, in particular since we have v > 0. In view of Lemma 

\'A.2\ we have for t S (0, v), 

A,{t + z) = A(t) + D(z), D.(t + z) = D(t) + D(z) 

this contradicts the fact that z < 00 and concludes the proof. □ 

Let £ C I])o(IK+ ) X Bo(M^) x B{R^) as defined at the beginning of Section and .F = 
B(Mf ) X ID)o(IRf )• 

For Net G i5 and (A, D) € .F, we define the function 

G(Net,A,D) = ||(A-r(D,Net),D-$(A,Net))||. 

The function G is continuous and such that 



G(Net, A,D) = 0^ 




= r(D,Net), 
= $(A,Net). 



Let PNet be the subspace of £ of piecewise linear networks: namely Net = (S, P, N) G I'Net 
if the functions u 1— > 'N^^\u),u S^^\u) and u 1— > P^*'-'^(n) are piecewise linear non-decreasing 
functions such that p{P{t)) < 1 for alU > and N^^) = for i ^ 5. We denote Net = (S, P, N). 

We proved that 

VNet G pNet, G(Net, A,D) = 4^ (A,D) = ^'(Net), 

where ^ has been explicitly defined above. We are exactly in the framework of Section^ In 
the next section we construct the mapping 5 : £" x — > I^^et • 



large deviations for queueing networks 
4. Sample path large deviations 
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In order to simplify the notations, we assume that N„(0) = for all n. This condition can 
be weakened to the standard condition: 

lim ilogP(N„(0) > e) = 0, 

n— >oo n 

for all e > 0. In this case, we have \^{x) = oo for all x 7^ and l''(0) = 0. 

It is possible to deal with the case where the initial condition satisfies a LDP as assumed in 
Theorem I2.1l bv using a standard conditioning argument (as done in !2T1 for example). 

4.1. Construction of the approximating sequence. This section is devoted to the proof of 
the following proposition: 

Proposition 4.1. We consider Net = (S,P,N) G S such that I^''*(Net) < 00 and such that 
there exists (A, D) G that satisfies the fixed point equation Ii0.1\) given by, 

A = r(D,Net), 
D = $(A,Net). 

There exists a sequence {Net„,}„ = 5(Net, A,D) such that 

(4.1) Net„ G PNet for all n; 

(4.2) Net„ Net; 

(4.3) ^'(Net„) ^ (A,D); 

(4.4) l^^*(Net„) ^ l^^*(Net). 

First note that since l^''*(Net) < 00, each process S, P, N is absolutely continuous and Net is 
well-defined. Moreover thanks to Lemma 15.31 the processes A and D are absolutely continuous 
too. 

The idea to construct the sequence {Net„}„ is to consider the piecewise approximation of 
the fixed point equation (|0.1|) . First consider the routing equation A = r(D,Net) for times t 
such that G N, 

K 

A«(i + 1/n) - A«(t) = N«(t + 1/n) - N»(t) + V ^^J'^)(D(-'')(t+)) {T>^^\t + l/n) - T>^^\t)), 
" ' ^ ^ ' 7^1 ^ ^ ' 

Ai''(A)(t) AW(N)(t) A(/'(D)(t) 

where we define the piece-wise linear process P^ ' ''{t) as follows, for s G (D(j)(t), D(j)(t + 1/n)), 

^o-,Or ^ = pO^-)(D(^-)(t + l/n))-PO^-)(DO-)(t)) 
" D(i)(t + l/n) -DO-)(t) 

if D(-'')(t + 1/n) / D(j)(t), and we take Pn '*^(D^-'H^)) = otherwise. In other words, we have 

p(f'^)(D(j')(t + 1/n)) - p(,^'^)(D(^')(t)) = ^(f'^)(D(^')(t+))(D(^')(t + 1/n) - T>^^\t)) 

= p(^'*)(D(j')(t + 1/n)) - p(j'*)(D(^')(t)) 

Note that {Pn'^\t)}i^j G since we have by the definition of 

^p(i,i)(DO-)(t + 1/n)) - p(j'^)(D(-'')(t)) < Ti^^\t + 1/n) - Ti^^\t), 
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but the matrix (Pn'^\Y)^^\t+)))ij may not be of spectral radius less than 1. 

To circumvent this difficulty, we modify slightly the processes as follows, (the variables r],en,S 
will be made precise latter) 

(4.5) aW(A) + ^ = aW(N) + — 

n n 

K 



3=1 \ 



n 



where we omit to write the time t and use the simplified notation Pn'*^ = Pn'*^(D(-'^(t+)). 

We have to find rj,en,6 such that (|4.5|) holds with rj^'^\ en\ 6^^^ non-negative and (5*-*^ = for 
i ^ S. These constraints are satisfied by the following choice: first take 6 such that J*-*^ > for 
all i G 5 and (5^*^ = for i ^ S. Let ri{5) = rj he the unique solution in M;^ of the following 
equation (recall that p{R) < 1), 

K 

j=i 

Note that 77*^*^ > for all i thanks to (|2.4)) . Finally let define e„(5) = as follows e 



n 



^ — - G (0, 1] (note that ei*^ = 1 if and only if A^^^(D) = 0). 

(D)+?7(») 

It is easy to see that (|4.5() holds since we have 

(1 - 6(f)) ( A(f)(D) + ^]= A(^-)(D), or, e^f) ( Ali^B) + U !^!^ 



which imply respectively that 



A»(A) = A»(N)+J^(l-6(^-))p(^'*) (^A(f-)(D) + ^) and. 



and summing these two equalities gives (|4.5j) . 

For 5 fixed, we define for s G (D(j)(t) + tr/((5), D(j) (t + 1/n) + (t + l/n)T/((5)), 

Pn? (^) = (1 - 6(f))^(^'^)(D(^-)(t+)) + 

where 6^(5) is defined as above. In view of Lemma l5.41 the matrix P^'^''(s) is of spectral radius 
less than one since e^ '' > for all j. Then as a direct consequence of 1)4. 5p . we have for G N, 

K 

(4.6) A»(t) + tr?(5) = N»(t) + t5 + Y,V^if {T>^^\t) + try(,5)). 
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If N„^5 is the polygonal approximation of t ^ + t6 with step 1/n, we have clearly 

Nn,5 ^ N + 5 as n tends to infinity. Similarly, we have as n tends to infinity, 



' (l-eW(t))pO-,0(DO-)(t))(Dl^toM 

+e(i)(t)i?(i,i)(D(i)(t) + ^(5)) if D(i)(t) > 0, 

R0''^)ri{6) otherwise. 



where e^^\t) = rj^^\6)/{rj^^^6) + ti^^^t)) < 1. Hence when n tends to infinity and 6 tends to 
zero, we have Nn,5 N and P^^^f P(J'*). 

We consider now the queueing equation D = ^{A, S) and construct the approximating se- 
quence for S. 

We begin with a first general lemma: given three processes A < D and S, we construct a 
piecewise linear function S„ (with step 1/n) as follows (with nt € N): 

• if A(t) = D(t) and A(t + l/n) = D(t + l/n), then S„(t + l/n)-S„(t) = S(t + l/n)-S(t); 

• otherwise, S„(t + 1/n) - S„(t) = D(t + 1/n) - D(t). 

We will denote this construction by S„ = T„(A,D, S). 

Lemma 4.1. Let (A, D, S) be absolutely continuous functions o/C(M:^) x Co(M:^) x Co{^^)such 
that <1>(A,S) = D. We denote Sn = T„(A,D,S). We have D„ = $(A„,S„) where (A„,D„) is 
the polygonal approximation of (A,D) with step 1/n and we have the following convergence as 
n tends to infinity: S„ — > S, S„ ^ S and \^ {Sn{t))dt \^{S{t))dt . 



Proof. We denote D„ = $(A„, S„). From the proof of Lemma ESI we have 

D„(t + 1/n) - D„(t) = inf |s„(t + 1/n) - S„(s) + A„(s) - D„(i)| A {S„(t + 1/n) - S„(t)} , 

t<s<t+l/n J 

since all the functions are linear on the interval (t, i + 1/n), we have (with nt G N), 

D„(t + l/n) -D„(t) = {A„(t + l/n) -D„(t)} A{S„(t + l/n) -S„(t)}. 

If t)n{t) = E)„(t), then we have clearly D^(t + 1/n) = D„(t + 1/n) since 

• if An{t) = Unit) and A„(t + 1/n) = D„(t + 1/n), then we have S{t + 1/n) - S{t) > 
D„(t + 1/n) - D„(t) = An{t + 1/n) - D„(t) see ((5TT|) for the inequality; 

• otherwise, S„(t + 1/n) — S„(t) = D„(t + 1/n) — D„(t) by definition and An{t + 1/n) > 
Bnit + l/n). 

This proves the first part of the lemma. Moreover it follows directly form the definition of T 
that S„ {t + ^) — Sn(t) < S (t + ^) — S{t), hence we have for all t, limsup^^o^^ S„(t) < S{t) by 
a continuity argument. The fact that S.„ S follows directly from Fatou's Lemma and the fact 
that Sn S. We now prove this last fact, let C = {t, A{t) = D(t)}. C is a closed set and 
according to Lemma l5.3| we have for all t € (the complementary set of C), S{t) = D(i). For 
such t G C"^, we have for e > sufficiently small and for sufficiently large n, An{u) ^ D„(n) for 
all |u — t| < e. Hence we have S„(t) = D„(t) D(t). Now for t € C° in the interior of C, we 
have clearly S„(t) -^S{t). Hence we have S„(t) S{t) for t e C" U C. 
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We prove the last statement of the lemma. Since any open set of R is a countable union of 
disjoint intervals, 

/ \^{Snit))dt= \^{'Dnit))dt < / \^{'D{t))dt, by Jensen's inequality 
Jc Jc Jc 

\^{S{t))dt, 



and also directly still by Jensen's inequality f' {Sn{t))dt < J^^ \^{S{t))dt. The convergence 
then follows from 

liminf / \^{Sn{t))dt > / liminf l^(S„(t))(it 

/•OO 

> / \^{S{t))dt, 
Jo 

where the first inequality is due to Fatou's Lemma and the second one to the lower semicontinuity 
ofl^. □ 

We define the sequence Net„ ,5 = (S„^5, P„^5, N„ ,5) where Sn,5{t) = T„(A(t) +r/t, D(t) + 
r]t, S{t) + r]t). Note that we have D(t) + r/t = $(A(t) + rjt, S{t) + rjt), hence Lemma UTTl applies . 
in particular, we have Sn,5{t) S(t) + r]{S) as n tends to infinity. 

We have Net„ ,5 G I'Net by construction and the sequence {Net„ 5,^}^ satisfies ()4.2() for some 
6n — 0. Moreover, we have thanks to (|4.6() and Lemma l4.ll 

J ^n,5 — r(Dn A, Net„ a), ^ / a t-» ^ iTr/TVT„j- ^ 
|d„:, = ci,(A„;,,Net„;,), ^(An,.,D„,,) = ^(NeW), 

where A„ ,5 and D„ ,5 are the polygonal approximation of A{t) + rjt and D(t) + r]t with step 1/n 
and ^ has been defined in Section |21 

For n ^ 00 and 5 — > 0, we have (A^^^jD^^^) — > (A,D), hence we have ^'(Net„.5) — > (A,D), 
i.e. the sequence {Net„_5^}„ satisfies (|4.3|) . 

We now show that (|4.4() is also satisfied. We fix T > and prove first that we have, for 6 
sufficiently small, 

(4.7) r l^(S(t)) + D{Pit)\\R) + |N(N(t))dt - er{5)T 
Jo 

< liminf / l^(S„,5(t)) + D(P„,5(t)||i?) + |N(N„,5(t))dt 

< limsup r l^(S„,5(t)) + D(P„,5(t)p) + \^{N^4t))dt 

n^oo Jo 

(4.8) < r \^iSit)) + DiP{t)\\R) + \^(Nit))dt + eri6)T, 

Jo 

where er{6) tends to zero as 5 tends to zero, from which 1)4. 4|) follows by monotonicity. 

We first deal with the case of the sequence of processes {Sn,s}n (we can restrict ourselves to 
the one dimensional case). We denote Ss{t) = S(t) + r]{5)t. 

We define ? = esssup{S(t), t < T} = inf{u, Leb[t < T, S(t) > u] = 0}, where Leb is for the 
Lebesgue measure. Since \^{S{t))dt < 00, <; belongs to the domain of 1^ which is open on the 
right. Hence we can find e > such that ? + e still belongs to this domain and take 5 such that 
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r]{S) < e. Moreover, since 1^ is convex, it is uniformly continuous on [0,<^ + e]. Hence, we can 
assume that we have (3{a) — s- as a ^ such that, 

yx,y£ [0,? + e], \x-y\<a^ \\^ {x) - \^ {y)\ </?(«). 
From Lemma l4.H we have 

\^{Sit))dt-Pir]{5))T < 



hm r\^{Sn4t))dt = r\^iSsit))dt 

"-*<^Jo Jo 



(4.9) < r \^iS{t))dt + Pir^i6))T. 







s 



Hence we proved (|4.7|) and H4.8|) for I 

The same kind of arguments can be repeated for ^n,5 which is just the polygonal approxima- 
tion of t N(t) + 5t. Note that {Net„,5}„ G V^^^ imphes that 'N^^^si*) = for all i ^ S. For 

i (z S, we can use the fact that the domain of 1^^*' is open as previously. In the case of Pn,(5, 
we can not use the argument on the openness of the domain, but we have D(i?W||i?W) = and 
then the convexity of D directly implies that D(P^*'*^||i?(*'') < □(P^*-' Hi?*^*)), from which we derive 
an equivalent of (|4.9|) . 

4.2. Exponential tightness. We first recall some definitions. A sequence of random variables 
{Xn}n G (M-^)^ is exponentially tight if 

lim limsup i logP(||X„|| > M) = -oo. 

For 5 > and T > 0, define the modulus of continuity in I}{E) by 

i(;'(X, T) := inf max sup d(X(s), X(t)), 

* s,te[ti-i,ti) 

where the infimum is over {ti} satisfying 

= to < ti < • • • < tm-i <T <tm 

and mini<j<„(tj - > 6. 

Theorem 4.1 of tells us: let Tq be a dense subset of M+. Suppose that for each t E Tq, 
{X„(t)}„ is exponentially tight. Then {X„}„ is exponentially tight in D(ii^) if and only if for 
each e > and T > 0, 

(4.10) limlimsup-logP('u;'(X„,5,r) > e) = -oo. 

<5^0 n^oo n 

A sequence of stochastic processes {X„}„ that is exponentially tight in 0{E) is C-exponentially 
tight if for each 77 > and T > 0, 

(4.11) limsup-logP(supd(X„(s),X„(s-)) > t]) = -00. 

Then Theorem 4.13 of ^] gives: an exponentially tight sequence {X„}„ in 0{E) is C- 
exponentially tight if and only if each rate function I that gives the LDP for a subsequence 
{'^n{k)}n{k)j Satisfies l(x) = 00 for each x € D(ii^) such that x ^ C{E). 
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The stochastic assumptions of Section 12.21 ensure that the sequence of processes {Net„}„ 
satisfies a LDP with good rate function (this imphes that the sequence is exponentially tight) 
giving an infinite mass to discontinuous path. Hence the sequence of processes {Net„}„ is 
C-exponentially tight. 

We have to show that the sequence of processes {(A„,D„)}„ is exponentially tight. The fact 
of dealing with non-decreasing processes simplifies the definitions. For X G D(Mf ) (or B(M^)) 
non-decreasing, 5 > and T > 0, we define ws(X.,T) = sup^gp^p] ll^(^ + 5) — X(t)||. We have 
clearly w'{'X.,6,T) = ws(X.,T) and if {X„(0)}n is exponentially tight then ()4.1U() implies that 
{X„(t)}„ is exponentially tight for each t > 0. Lemmas 15. II and 15.21 show that conditions (|4.1U() 
and (|4.11j) are satisfied for the sequence of processes {(A„,D„)}„. The exponential tightness of 
{(A„(0),D„(0))}„ is clear since A„(0) = D„(0) = 0. 

4.3. Large deviations results. 

Proposition 4.1. The sequence of processes {(A„,D„)}„ satisfies a LDP in D(M:^) x D(R:^) 
with good rate function \^'^ . For A,D absolutely continuous and such that A(0) = D(0) = 
and A > D, 1"^'^ is given by 

/•oo 

(4.12) l^'°(A,D)=/ //(A(s),D(s),A(s),D(s))ds, 

Jo 

where H{A, D, A, D) := infp^Ar h{A, D, A, D, P, N), with h given by, 
h{A,D,A,D,P,N) := 

E i'''(^^^^)i0«>,«}+ E i^"^(i>«) + E^^^^D(p»p») + |N(iv) 

ieE{A,D) i^E{A,D) i 

where E{A, D) = {i, A^*) = L'W} and with the infimum taken over the set of {P, N) G M-^ x M.^ 
such that 

A = N + P^b. 
For all other A, T), we have |-'^'^(A,D) = oo. 

Proof. Thanks to the results of previous sections, conditions of Proposition 11.11 are satisfied and 
we define 

(4.13) P'D(A,D) = inf I hm |N-t(Net„), {Net„}„ G 5(A,D)| , 

where we recall that 5(A,D) = UNet'5(Net, A, D), and 5(Net, A,D) is defined in Proposition 
Ol We have to show that T^'° = given by ^J^. 

Consider Net G I^Net and let (A,D) = ^(Net). Let r = {0 = to < *i < • • • } be such that 
the processes A,D,S,N and D o P have a constant derivative on each (tk,tk+i). Then from 
A = r(D,Net), we derive 

AW(t) = N^'\t) + ^-b^^\t)P'-j''\-D^^\t)). 
j 

From D = $(A, S), we get the following constraints: 

• if AW(tfc) > DW(tfc) or AW(tfe+i) > DW(tfe+i), then we have DW(t) = S^^{t) for 
t G {tk,tk+i)', 

• otherwise AW(t) = DW(t) for t G (tfc,ifc+i) and we have SW(t) > AW(t) = ti^'\t) for 
t G (tfc, tfc+i). 
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Now we can compute l^''*(Net) as follows 

i€E(A,D) i<^E{A,D) 

^D(^')(s)D(p(^')(s)||i?(^'))ds 
i 

> / /i(A(s),D(s),A(s),D(s),P(s),N(s))(is > l^'°(^(Net)), 

since for i G £;(A(s), D(s)), we have l^*''(S«(s)) > 1^'"' (D(^)(s))1|d(,)(^)>^(,)j because SW(s) > 

D^*^(s) and 1^''' is non-negative, convex with /i*-*-* as unique zero. Hence, we have I^'-^ > |-^'^. 

Consider now (A, D) such that l^'-^(A, D) < oo, then we denote by (p(s), n(s)) the argument 
that achieves the minimum in ii'(A(s), D(s), A(s), D(s)) for any fixed s (note that /i is a good 
rate function). Let P(D(t)) = p(s)(is and N(t) = n(s)ds, note that p and n are measurable 
since H \s a, good rate function. We have A = r(D, P, N). Now define s(s) as follows: 

• if A«(s) = dW(s) then s«(s) = D«(s) V 

• if AW(s) > DW(s) then sW(s) = DW(s). 

We have D = ^>(A,S) with S(t) = s(s)ds. Hence we have (A,D) = (r(D, Net), $(A, Net)) 
for Net = (S,P,N) and l'^^*(Net) = l'^'^(A,D) < oo by construction. Hence the sequence 
5(Net,A,D) = {Net„}„ is well-defined and we have T^'°(A,D) < lim^^oo l^'^HNetn) = 
|Ne*(Net) = l^'°(A,D). □ 

From this proposition, it is quite easy to derive a LDP for the process Qn(t) := A„(t) — D„(t) 
counting the number of customers in each queue. Thanks to the contraction principle, we have 

|Q(Q) = inf{l^'°(A,D),Q = A-D}, 

which gives directly Theorem 12.11 

Acknowledgements. I am thankful to Anatolii Puhalskii for insightful comments and for 
providing me with a copy of 

5. Appendix 

5.1. Properties of the map T and For X € ©(Mf ), 5 > and T > 0, we define Wi{^, T) = 
SUPtg[0,T] l|X(t + 5) -X(t)||. 

Lemma 5.1. We have 

ws{HX,Y),T) <ws{Y,T). 

Proof. It is clearly sufficient to consider the case K = 1. We will prove that 
(5.1) $(X,Y)(t + 5) -$(X,Y)(t) <Y{t + 5) - Y{t), 

from which the lemma follows. If $(X, Y)(t) = Y(t), then we have $(X, Y)(t + 5) < Y{t + 5) 
and (|5.1|) is clear. 

Assume now that $(X, Y)(t) = info<s<t {Y{t) - Y(s) + X(s)} < Y(t). We have 
Y{t + 6)- Y{s) + X(s) = Y{t) - Y{s) + X(s) + Y{t + 5)- Y{t), 
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and (|5.H) follows by taking the minimum in s E [0,t] and observing that $(X,Y)(t + 6) < 
info<.<t{Y(t + J) -Y(s) + X(s)}. □ 

The following lemma is clear: 

Lemma 5.2. We have 

wsinx,P,N),T) < ws{N,T)+ws{P, ||X(r)||). 

Lemma 5.3. Assume S S Bo(M+) is absolutely continuous, then for any A € ][!)(M+), we have 
D := <I>(A,S) is absolutely continuous and, 

• for all t such that A(t) > D(t), we have D(t) = S{t); 

• if A(t) = D(t) fort G {u,v) with u <v, then we have S{t) > A{t) = D(t) fort £ {u,v). 

Proof. It follows directly form (|5.1|) that if S is absolutely continuous, then <I>(X, S) is absolutely 
continuous for any X. The rest of the lemma is obvious. □ 

5.2. Auxiliary results. 

Lemma 5.4. Given a substochastic matrix R such that p{R) < 1 and a substochastic matrix P 
such that the support of P is included in the support of R, i.e. R^^'^^ = =^ p(*.i) = Q. Then for 
any e such that < e*^*) < 1 for all i, the matrix with coefficients M^*'-') = (1 — e(*))p{*'J)+e(*)i?(*'-?) 
is of spectral radius less than 1 . 

Proof. By a suitable permutation of rows and columns, we can assume that R is given in its 
canonical form 



(5.2) R 



( Si{R) * * * \ 

S2{R) * * 

■•. * 

\ SniR) J 



where each Si{R) is an irreducible matrix. We have p{R) < 1 if and only if each Si{R) is not a 
stochastic matrix. 

In view of the assumption on the support of P, the matrix P has the same structure as (|5.2|) 
and we have with the same notation as above, Si{M) which is an irreducible and not stochastic 
matrix. □ 

5.3. An example. In this section, we construct 2 different sequences of Jackson networks Net^ 
and Net^ such that their fluid limits are the same 

Net^ Net and Net^ Net, 

but such that 

(A^,D^) = M/(Neti) ^ (Ai,Di), 

(A2,D2) = M/(Net2) ^ (A2,d2), 

with (A\Di) / (A2,D2). 
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We consider a toy example with only one station (hence we omit the superscript .^^^ that 
refers to that only station). Once a customer is served, he can either go out of the network or 
go back to this same node. We define the following driving sequences: 

= (l^_^,n,l^_^,n,...), 

n n 

= a(l,l,...), 
with a < 1. We define now two different routing sequences 

= (2_^,1_J,...), 

n+l n+1 

\xn\ 71— \xn\ \_xn\ n~ \xn\ 

where x < 1. We denote by Net^ = {a", i/", T"} and Net^ = {a", i/"(x), T"}. z^"(x) is obtained 
from v'^ by only interchanging a 1 and a 2. Hence we have 

Net^ ^ Net and Net^ ^ Net. 

Indeed the fluid network Net is given on Figure H 




Figure 1. Fluid networks: Net 



In the fluid limit, in case 1, the queue is always empty and the departure process is the same as 
the arrival process N. In case 2, the fluid limit of the departure process and the queue length 
process is given on Figure [21 




xa 



Figure 2. Departure process and queue length process 
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To explain D^, we write for each arrival (number on the left) the couple corresponding to: 
the inter-arrival time | the routing decision (1 means that the customer goes back in the queue 
and 2 means that the customer leaves the network): 



1 
2 
3 

[xn\ 
[xn\ + 1 
[xn] + 2 

n 
n + 1 

n + 2 

n + 3 



1 
1 
1 

1 
1 
1 

1 

n 

1 
1 



2 
2 
2 

2 

1,2 
2 



1,...,1,2 
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